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ABSTRACT 

We present an updated catalogue of M31 globular clusters (GCs) based on images 
from the Wide Field CAMera (WFCAM) on the UK Infrared Telescope and from the 
Sloan Digital Sky Survey (SDSS). Our catalogue includes new, self-consistent ugriz and 
K-band photometry of these clusters. We discuss the difficulty of obtaining accurate 
photometry of clusters projected against M31 due to small scale background structure 
in the galaxy. We consider the effect of this on the accuracy of our photometry and 
provide realistic photometric error estimates. We investigate possible contamination in 
the current M31 GC catalogues using the excellent spatial resolution of these WFCAM 
images combined with the SDSS multicolour photometry. We identify a large popula- 
tion of clusters with very blue colours. Most of these have recently been proposed by 
other work as young clusters. We distinguish between these, and old clusters, in the 
final classifications. Our final catalogue includes 416 old clusters, 156 young clusters 
and 373 candidate clusters. We also investigate the structure of M31's old GCs using 
previously published King model fits to these WFCAM images. We demonstrate that 
the structure and colours of M31's old GC system are similar to those of the Milky 
Way. One GC (B383) is found to be significantly brighter in previous observations 
than observed here. We investigate all of the previous photometry of this GC and 
suggest that this variability appears to be genuine and short lived. We propose that 
the large increase in its luminosity my have been due to a classical nova in the GC at 
the time of the previous observations in 1989. 

Key words: galaxies: individual: M31 - galaxies: star clusters - globular clusters: 
general 



1 INTRODUCTION 

Globular clusters (GCs) are among the oldest known stel- 
lar systems. They typically have ages similar to those of 
their host galaxies, making them ideal probes into galaxy 
formation and evolution. The properties of GCs vary signif- 
icantly. However, individual clusters contain populations of 
stars with similar ages and metallicities. This makes them 
unique locations for studying stellar evolution. 

The Milky Way's GCs still represent the best studied 
GC system. While the study of these clusters has led to many 
advances, the Milky Way contains relatively few GCs (~150 
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GCs: lHarri3ll996l ). many of which have high foreground ex- 
tinction, making them hard to study. By determining the 
properties of extragalactic GCs, we are able to study a more 
diverse population and ensure our current conclusions are 
not biased by the Milky Way's clusters being atypical. 

For extragalactic GCs, it is very difficult to resolve in- 
dividual stars in the clusters. However, it is possible to esti- 
mate many properties of a GC from its integrated light. For 
example: the masses of GCs can be estimated by assuming 
a mass to light ratio; combined optical and near infrared 
colours of GCs can be used to (at least partially) break the 
age and metallicity degene r acy and estimate these parame- 
ters (e.g. IPuzia et al ] |2002l : Ijiang et al.l[2003l : iHempel et al.l 
and their structural paramete rs can be estimated 
by fitting their density profiles (e.g. iBarmbv et al.l l2007l : 
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Jordan et al] l2007l : iMcLaughlin et alj l2008l : iPeacock etal] 
20091 ). The colours of GCs and GC candidates are also very 



useful in selecting genuine GCs from stellar asterisms and 
background galaxies. Good multi-wavelength photometry of 
GCs is therefore highly desireable. 



1.1 The M31 GC system 

The proximity of M31, and its relatively lar ge GC popula- 
tion compared with the Milky Way (~400: iBarmbv et al.l 
l200d '). makes it the ideal location to study extragalactic 
globular clusters. Its clusters have been the focus of m any 
studies d a ting b ack to the early work of lHubbl^ ()l932h and 
IVetesnikI l| 19621 ). However, attempts to study its clusters 
have faced several challenges. Photometry of these clusters 
is complicated by many of them being projected against the 
bright and non-uniform structure of M31 itself. The galaxy's 
proximity also results in the GC system extending over a 
wide region of the sky, with cluste rs recently fo und beyond 
4° from the centre of the galaxy (|Huxor et al. 2008). This 
means that surveys with large fields of view are required in 
order to study the GC system. It is also difficult to confirm 
GCs in M31 based on spectroscopy alone. This is because 
the velocity distribution of its GC system overlaps that of 
Milky Way halo stars. 

Over the past decades there have been several 
large catalogues of M31's GCs i ncludi n g those of: 
iBattistirii et al.l lll987h: iBarmbv et al.l (|2000l ): lGalleti et all 
■ (|2 



l|2004h : iKim et al.r l|2007^ . In addition to these cata 



posed (e.g. iBattistini et al.l Il993l: iMocheiska et al 
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made considerable progress in removing contamination from 
the cluster catalogues due to either background galaxies 



(e.g. lRacinjri 991: Ba rmbv et alll200d: IPerrett et al.ll2002l : 



iGalleti et ai]r2004 : .Kim et al.ll2007l : ICaldweU et "al]l2009[ ) or 

stars an d asterisms both in the Milky Way and M31 it- 
sclf (e.g. Barmbv ct al. 2000:'Cohen. Matthews, fc CameronI 



12005 ; Huxor ct al. 2008; Ca ldwell et al . 2009). However, de- 
spite this work, it is likely that there remains significant con- 
tamination in the current catalogues of M31 clusters, espe- 
cially at the faint end of the GC luminosity function. These 
studies have also resulted in a large num ber of unconfirmed 
candidate clusters (currently over 1000: iGaUeti et al.ll2004l ) 
whose true nature remains uncertain. 

It has been known for many years that some of the 
proposed GCs in M31 have very blue colours. Recent work 
has identified that a large number of the c lusters in the 
current catalogues are young clusters (e.g . iBeaslev et al.l 
l2004l : lFusi Pecci et alllioosl : iRev et al.ll2007i ). A comprehen- 
sive catalogue of young clusters in M31 ha s recently been 
publis hed from the spectroscopic survey of ICaldwell et ahl 
lj2009l ). Compared with the young open clusters in the Milky 
Way, these clusters have relatively high masses (< 10^ Mq), 
akin to the young clusters observed in the Large Magel- 
lanic Cloud. A recent Hubble Space Telescope (HST) study 
of 23 of these young clusters suggested that on average 
they are larger and more c oncentrated than typical old 
cluster s (IBarmbv et al.l 20091 ) . Most of the clusters stud- 
ied by IBarmbv et al.l ( 20091 ) are found to have dissolution 
timescales of less than a few Gyr and are therefore not ex- 



pected to evolve into typical old globular clusters. Whether 
these clusters are massive open clusters, young globular clus- 
ters or a mix of both, it is clear that they represent a dif- 
ferent population to the classical old GCs also observed in 
M31 (which are the focus of this study). We therefore dis- 
tinguish between the two populations in our classifications 
and conclusions. 

While these previous studies have provided a wealth of 
information on the M31 GC system they have also resulted 
in a rather heterogeneous sample. For example, the excellent 
and commonly used R evised Bolog na Catalogue (hereafter 
RBC) of M31 GCs by Galleti e t al. (2004) includes photom- 
etry from many different authors using different telescopes 
and in some cases differen t (homogenised ) filter s. This has 
been previously noted by I Caldwell et al.l (|2009l ) (hereafter 
C09) who recently published new V-band photometry for 
many RBC sources which were locate d in the Local Grou p 
Galaxy Survey images of M31 (LGGS: iMassev et al.ll2006l ). 
While this work provides excellent deep V-band photometry, 
the survey does not cover the outer clusters and candidates 
and does not provide colour information. The most com- 
plete set of optical colours, derived in a consistent manne r, 
are still that of the Barmby catalogue (|Barmbv et al.ll2000l ) . 
This work presented self-consistent UBVRI colours for many 
of their clusters. However, it is incomplete in some of these 
bands and only provides new photometry 285 clusters. For 
these reasons we chose to produce new, self-consistent, op- 
tical photometry for the proposed GCs and GC candidates 
in the RBC using images from the Sloan Digital Sky Survey 
(SDSS). The excellent calibration and large field of view of 
this survey is ideal for studying such an extended system. 
Details of this photometry are presented in section 2.2. 

The study of M31's GCs in the near infrared (NIR) is 
very useful both for confirming genuin e GCs and for esti- 
mating th eir ages and mctallicit ies ( e.g. IBarmbv et al. I I2OOOI ; 
iGalleti ei al.l |2004| : [Fan et all 12009 ). The first maj o r sur - 
vey of M31's GCs in the NIR was bv IBarmbv et al.l (I2OO0I ) 
who used pointed observations of individual clusters to ob- 
tain K-band photometry of 228 of their clusters. More re- 
cently |GalietLe£^L| (|2004h obtained NIR photometry in the 
J,H and K-bands of 279 of their confirmed GCs from the 
2 Micron All Sky Survey (2MASS). The spatial coverage 
of 2MASS makes it ideal for such a project. However, the 
survey is relatively shallow and has relatively poor spatial 
resolution. We have obtained new deep K-band photometry 
using the Wide Field Camera on the UK Infrared Telescope 
to determine the K-band magnitude of M31's GCs across the 
entire GC luminosity f unction. Some results of this survey 
are already published in lPeacock et al. I l|2009l ). In addition to 
providing the first K-band photometry for 126 GCs marked 
as confirmed in the RBC, the excellent spatial resolution 
of these images is very useful for removing stellar sources 
from genuine clusters, and for investigating the density pro- 
files of the clusters. Details of this new K-band photometry 
are presented in section 2.3, while the classifications of the 
proposed clusters and candidates are considered in section 
3. 

The proximity of M31 makes it the ideal location for 
studying the structure of extragalactic GCs. Determination 
of the size and density of GCs is very useful in investigat- 
ing: stellar evolution; galaxy formation and evolution; con- 
straining N-body simulations; and investigating exotic ob- 
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Figure 1. Coverage of the SDSS and WFCAM images used. For reference all objects listed in the RBC are shown in red. Only images 
covering the locations of these objects were extracted from the SDSS archive. The green ellipse indicates the D25 ellipse of M31. The 
grid represents 2° X 2° squares on the sky (13.6 X 13.6 kpc at the distance of M31) and highlights the spatial extent of the GC system. 



jects in GCs like X-ray binaries, blue stragglers and hori- 
zontal branch/extreme horizontal branch stars. The struc- 
tural parameters for some of M31's GCs have been measured 
using Hubble Space Telescope (HST): Faint Object C amera 
(FOC) images of 13 clusters l|Fusi Pecci et al.Hl994l '): Wide 
Field Planetary Camera (WFPC2 ) images of 50 clusters 
l|Barmbv. Holland, fc Huchral 120021 '): and Advanced Cam- 
era for Surveys (ACS) and Space Telescope Imaging Spec- 
trograph (STIS) im ages of 15 and 19 cluster s respectively 
l|Barmbv et al.ll2007l '). In lPeacock et al.l l|2009l ) w e presented 
the results of fitting PSF convolved King models (lKingj|l966l ) 
to the ground based WFCAM images used here in order to 
estimate the structural parameters for 239 clusters. In sec- 
tion 4 we discuss the structure of M31's GCs. 

2 PHOTOMETRY OF GCS AND CANDIDATES 
2.1 Identification of GCs 

In the following analysis we consider all the GCs and 
GC candidates listed the RBC (their class 1/8 and 2 
objects respecti vely). Based on the original catalogue of 
[Battistini et al] (Il987l ). this catalogue has been regularly 
updated to include the results from most new studies. This 
version of the catalogue (v3.5) include s the newly discov- 
ered GCs in the o uter regions on M31 (|Mackev et al.ll2006l : 
Huxor ot al. 2008) and the new GCs and candidates from 
Kim ot a l. (2007) (hereafter K07: their class A and B/C ob- 
jects respectively). We also consider the catalogue of C09 
which includes some additional clusters and gives updated 
locations and classifications for many of the objects in the 
RBC based on images from the LGGS or Digital Sky Survey 



and/or Hectospec spectroscopy. This combined catalogue is 
used to identify the known GCs and candidates in the fol- 
lowing analysis. 



2.2 Optical photometry 

2.2.1 ugriz data 

To obtain self-consistent optical photometry of M31's clus- 
ters and candidates we extracted images of M31 from the 
SDSS archive. Since M31 is at a relatively low Galactic 
latitude of -21°, it is not included in the standard survey 
field. However drift scan i mages of M31 were obtained by 
the SDSS 2.5m telescope (jAdelman- McCarthy et all 120071 ) 
in 2002 as part of a special run during a period when th e 
survey's primary field was not available (|Zucker et al.ll20o3 ). 
The runs used (3366, 3367, 6426 and 7210) provide image s 
in the five SDSS bandpasses {ugriz: iFukugita et al.|[l996l ). 
Each of the observations takes images in these bands simul- 
taneously meaning that they are taken under the same at- 
mospheric conditions. The seeing for different observations 
varied significantly between 1.1-2.1" in g (meaning that faint 
GCs could appear as point sources in some of these im- 
ages). The 3(7 detection limits of these images were verified 
to be similar to the standard survey (u < 22.0, g < 22.2, r < 
22.2, i < 21.3, z < 20.5). These data were found to cover and 
detect 92% (gri), 90% {z) and 73% (u) of the 1558 clusters 
and candidates in the current RBC. Two of these GCs were 
saturated in the r and i bands and one was saturated in 
the (?-band. We do not provide new photometry for these 
clusters but good photometry is already available for these 
very bright clusters from previous studies. 
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Figure 2. Difference between tlie location of objects in our im- 
ages and their location in the RBC (solid) and C09 (open). In 
both cases, data are grouped into 0.2" bins. 



the SDSS images may not appear extended due to the poor 
angular resolution of some of the images. Also, we wish to 
provide photometry for potentially misclassified stars in ad- 
dition to the extended clusters. The M31 GC catalogue was 
therefore matched to sources with stellar profiles in addition 
to those with extended profiles. In the few cases where mul- 
tiple sources were located within 3" of the quoted location, 
priority was given first to sources flagged as extended and 
then to the closest source to the quoted location. 

Figure [2] shows the difference between our positions and 
the positions quoted in RBC (solid) and C09 (open). We find 
excellent agreement between our locations and those of C09 
but find that the difference in the positions of many sources 
in the RBC are greater than 1". The errors in the positions 
of some sources in the RBC were noted and discussed by 
C09. We note the strong agreement between our locations 
and those of C09 and use their locations to identify GCs and 
candidates. 



We extracted all images covering the locations of con- 
firmed and candidate clusters from the SDSS Supplemen- 
tal Archive. Figure [T] shows the coverage of these data 
and demonstrates that most known clusters and candidates 
(red circles) are covered. These imag es have been processed 
through the standard SDSS pipeline (jStoughton et al]|2002l ) 
which both reduces the raw images and produces a catalogue 
of sources in each image. Since the SDSS extraction and pho- 
tometry routines are not designed to work in crowded fields 
(like M31), the default catalogues produced by the pipeline 
can not be used for photometry of the clusters. Instead we 
performed photometry on the images as described in the 
next section. 

The photometric zero points for these images were 
calculated using the calibration coefficients produced by 
the pipeline. These calibratio ns place the magni tudes on 
the AB photometric system [Oke fc GunnI (|l983l ): the u- 
band zeropoint has previously been found to be slightly 
offset from the AB system by mab = msdss — 0.04 mag 
l|Bohlin. Dickinson. &: Calzetti|[2001^ ■ this correction is not 
applied to our photometry]. This calibration is known to 
give magnitudes accurate to ~0.01 mag. 

2.2.2 Identification and locations of GCs 

Catalogues of all sources in each of the uqriz images were 
prod uced using the program SExtractor l|Bertin fc Arnouta 
Il996t ). This detected and located every source in each fil- 
ter, performed initial aperture photometry, and gave an es- 
timate of the stellarity of each source based on the PSF 
of its host image. Sources were identified using a mini- 
mum detection area (DETECT_MINAREA) of 3 pixels and 
a detection (DETECT.THRESH) and analysis threshold 
(ANALYSIS.THRESH) of 1.5o-. This threshold was chosen 
to ensure that the majority of the clusters profiles were in- 
cluded in the analysis and to help to separate genuinely ex- 
tended objects from merged sources. 

This source catalogue was matched to our combined 
catalogue of known M31 GCs and candidates (described in 
section 2.1) based on astrometry. We identified all objects 
within 3" of the locations quoted in the RBC and (sepa- 
rately) to their locations in C09. Some genuine clusters in 



2.2.3 Photometry 

Photometry of all clusters and candidates was obtained us- 
ing SExtractor's simple aperture photometry. We also con- 
sidered using the IRAF:APPHOT routines to perform the 
aperture photometry but SExtractor was found to deal bet- 
ter with contamination from neighbouring sources. This is 
a significant problem when using aperture photometry to 
obtain magnitudes of extended sources in a crowded region 
like M31. To minimise the effects of neighbouring sources 
within the GC aperture, SExtractor masks all other sources 
detected in the aperture and replaces them with pixels from 
symmetrically opposite the source. 

For background estimation we considered the use of 
both local and global solutions. To produce a global esti- 
mate of the background SExtractor produces a smoothed 
background map for each image. We chose to create this with 
a BACK_FILTERSIZE of 3 and a BACK_SIZE of 64 pixels. 
By examination of the background maps produced by SEx- 
tractor, this method was found to give a good estimation 
of genuine background variation (due mainly to structure in 
M31 itself) without subtracting fiux from the sources of in- 
terest. This was compared with the photometry produced 
using local backgrounds (calculated around the isophotal 
limits of the sources). In most cases good agreement was 
found between the two methods. However the local back- 
ground estimates were found to deal better with the most 
strongly varying background regions (near the center of the 
galaxy and its spiral arms) . For this reason local background 
estimation was used for the final photometry. 

To determine the total luminosity of each cluster, we 
produced curves of growth from i^-band photometry ob- 
tained through apertures with a radii in the range 2.8 — 10.6" 
with 0.6" increments. These were used to determine the 
aperture size required to enclose the total cluster light. The 
best aperture was determined independently for each object. 
This method ensures that we measure the total cluster lu- 
minosity correctly for the largest clusters. While the use of 
smaller apertures for smaller clusters maximises the signal to 
noise and minimises the contamination from nearby sources. 
The aperture size used to determine the total magnitude of 
each cluster is quoted in table 1. The average aperture radius 
used was ^-^5.8", with 87% of the apertures ^8.2". The ugriz 
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colours of the clusters were measured through 4" apertures. 
We also measured the colours using the aperture determined 
for the total g-band magnitude. This confirmed that there 
were no significant aperture effects due to the use of smaller 
apertures. For the final colours we chose to use the smaller 
aperture size in order to maximise the signal to noise and 
minimise the contamination from nearby sources. 

The statistical errors across the GC luminosity func- 
tion are in most cases less than 0.05 mag in gn. In general 
the u and z bands have slightly larger errors as they have 
slightly lower signal to noise. However, there are additional 
systematic errors which need to be considered. 

Firstly, the errors on the zero point calculation are esti- 
mated at 0.01 mag. This error dominates over the statistical 
errors for many of the bright clusters. The other significant 
source of error for some of the clusters is due to contamina- 
tion from nearby sources and the error on the background es- 
timation. In the bluer wavelengths there is significant small 
scale structure in M31. For clusters projected against the 
densest regions of the galaxy, this makes background sub- 
traction difficult as it can vary on scales smaller than the 
cluster of interest. The issue of background estimation is 
found to be particularly significant for the g-band photom- 
etry. In this filter there is significant small scale structure 
across M31 (this structure is less significant in the shallower 
M-band images where the statistical errors are larger) . In or- 
der to estimate the error on our background estimation we 
repeated our photometry with apertures 1.8" larger than the 
aperture used for the total magnitudes. With perfect back- 
ground estimation, the determined fiux should be the same 
through both apertures (within the photometric errors) . The 
difference in luminosity of the cluster through each aperture 
can therefore be used to give an estimation of the error on 
the background estimation. In most cases this estimated er- 
ror is quite small, with a median value ~0.015 mag, but for 
a few clusters it can reach '^O.l mag. This additional error 
is combined with the calibration and statistical error and 
included in table 1 as <Tg,tot- 

The error on the ugnz colours should be less affected by 
these effects. This is because they are often taken through 
smaller apertures, also we expect that possible errors on 
the background level in each filter should, at least partially, 
cancel. For this reason the errors on the colours quoted in 
table 1 are only the statistical errors. 

2.2.4 Comparison with previous photometry 

To date the best set of optical colours of M31 GCs derived 
in a consistent manner is that of the Barmby catalogue. 
This includes colours for 285 clusters in the Johnson UB- 
VRI bands obtained through 8" apertures. The catalogue 
also contains photometry for an additional 160 clusters 
collated from other studies. This collated photometry is 



taken mainly f rom the work of : iReed. Harris, fc Harria 
iBattistini et al.l (119931): iMocheiska et al l 




Sharpy fc Liutv 



(11983): 



Sharov, Lyutyj, 

ISharov, Lvutvi. fc Esipov, (Ji}9 



fc Esipovl 
(199^; 



Sharov. Lvutvi . 
Ikonnikoval 

I). For fuU details of 
the sources and reliability of this addi tional photometry we 
refer the reader to the description in iBarmbv et al.l (|2000l ) 
and the references therein. 

There is little previous optical photometry in the ugriz 
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Figure 3. Top: comparison between our total cluster magnitudes 
and those from C09. Bottom: comparison between our total clus- 
ter magnitudes and colours and those from the Barmby catalogue 
(B2000). In all cases the y axis is our photometry minus that 
obtained previously. The highlighted point indicates the cluster 
B383. 



bands with which to compare our results. However it is pos- 
sible to compare our colours with those of the Barmby cata- 
logue by transforming between the UBVRI and ugriz bands. 
This was done u sing the following the transformations from 
Ijester et all ((2005) 



V ^g~ 0.59 x{g-r)- 0.01 ± 0.01 
u~g = 1.28 X (U - B) -I- 1.13 ± 0.06 
g-r = 1.02 X (B - V) - 0.22 ± 0.04 
r - i = 0.91 X (Rc - Ic) - 0.20 ± 0.03 



(1) 
(2) 
(3) 
(4) 



These trans f ormat ions are based on all stars studied by 
Ijester et al.l (j2005l ). Applying these to the colours of glob- 
ular clusters may introduce a slightly larger error than the 
quoted rms residuals as globular clusters are stellar popula- 
tions rather than single stars. However, they can be used to 
check for consistency with this previous work. 

Figure [3] compares our colours for confirmed clusters 
with the transformed colours from the Barmby catalogue. 
The errors quoted include the residual from the transforma- 
tions and the errors on our photometry only. The scatter is 
therefore expected to be larger than la due to errors on the 
previous photometry. It can be seen that reasonable agree- 
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Table 2. Photometry of B383 



Source of photometry 


Observation date 


Detector 


AU 


AB 


AV 


AR 


SDSS (This study) 


2002 October 06 


CCD 


[ 17.15 


16.70 


15.72 


15.06 ]* 


Sharov, Lvutvi, & Ikonnikova fl992) 


1990 October 15-18 


Photoelectric 


-0.13 


-0.06 


-0.06 




Reed. Harris. & Harris (1992) 


1989 August 23-30 


CCD 




0.54 


0.39 


0.18 


Sharov & Liutvi figSSl 


1980 October 8-13 


Photoelectric 


-0.05 


-0.03 


-0.06 




Battistini et al. f 1987) 


1977-1981 


Plate 


0.16 


0.07 


0.16 





The magnitude difference between the SDSS photometry presented here and that found by previous work. 
* The total magnitude of the cluster in the SDSS images (transformed to Johnson filter system) 



ment is found between the u-g and r-i colours of the clusters. 
For the g-r colours a slight offset of 0.035 mags is found. 
However, this is within the rms scatter of the transforma- 
tions. We believe this offset may be due to the errors in the 
transformations (due to the difference in the spectrum of a 
typical globular cluster compared with a single star) , rather 
than a genuine offset between the colours. We therefore be- 
lieve that for most clusters our colours are consistent with 
the previous UBVRI colours in the Barmby catalogue. 

The top panels of figure [3] compares the total magni- 
tudes of the clusters obtained here with V-band photometry 
from the Barmby catalogue and the more recent photome- 
try of C09. To compare the total magnitudes of the clusters, 
our (j-band photometry was transformed to the V-band us- 
ing equation [T] For most clusters good agreement is found 
between our magnitudes and those in the Barmby catalogue 
(again errors on the Barmby catalogue photometry are not 
included). However, it can be seen that there are some signif- 
icant outliers. Some of the brightest clusters are brighter in 
our photometry than found previously. This is likely due to 
our use of larger apertures for larger clusters. Many of these 
bright clusters are found to extend beyond the 8" apertures 
used to obtain the Barmby catalogue photometry. For the 
fainter clusters we identify a group of 7 clusters which are 
fainter than expected. These clusters all have nearby neigh- 
bours, the effects of which we attempt to remove from our 
photometry but believe are included in the previous pho- 
tometry. We therefore believe our values for these clusters 
to represent the actual cluster magnitudes better. Another 
group of faint clusters are found which are brighter than ex- 
pected. The majority of these clusters are in dense regions 
near the galaxy center or spiral arms and we believe the dif- 
ferences are due to errors in the background estimation. It 
is very difficult to estimate the background accurately for 
regions with variations on the scales of the clusters them- 
selves. It is unclear which photometry is more accurate for 
these few clusters, although our use of smaller apertures for 
smaller clusters should minimise this effect. In the Barmby 
catalogue they subtract light from the bulge of M31 before 
performing photometry using a ring median filter. We re- 
peated our photometry using a similar method but did not 
find significant differences in our photometry. Background 
estimation for clusters in these dense regions is an inher- 
ent problem in finding their absolute magnitudes. It should 
be noted that, while we attempt to account for this in the 
quoted errors on our photometry, the errors for some clusters 
in these dense regions may be larger than quoted. 

It can also be seen that excellent agreement is found be- 
tween between our photometry and that of C09. The errors 
on this comparison are larger due to the inclusion of the er- 



rors quoted by C09. The group of clusters which were fainter 
in our photometry than the Barmby catalogue are found 
to agree well with this photometry. This is likely due to 
C09 also subtracting the effects of nearby sources from their 
photometry. They also use a similar method of increasing 
their aperture size for larger clusters, and our photometry 
for brighter clusters agrees with theirs. We again identify a 
few clusters in dense background regions whose magnitudes 
are slightly fainter than expected. However, not all of these 
are the same outliers as those found in the Barmby cata- 
logue. This highlights the difficulty of accurately obtaining 
integrated magnitudes for clusters in these regions. 



2.2.5 Variability in B383: a classical nova? 

Figure |3] identifies one relatively bright cluster (B383) which 
is significantly fainter (AV = 0.39) in our photometry than 
found in previous photom etry. This cluster was not observed 
by Barmby et al.l (|2000l ) and its BVR band photometry in 
b oth the Barmby catalogue and the RBC are from the work 
of lReed. Harris, fc Harrij (|l992l ). This cluster has high sig- 
nal to noise, a relatively clean background, and its mag- 
nitude was obtained through a similar sized aperture to 
that used previously (7.8"). The cluster is present in two 
different SDSS observations and the magnitudes obtained 
from each agree very well. Table [2] compares our photom- 
etry with other previous observations of B383. It can be 
seen that there is good agreemen t between our photome - 



try and the previous photomet r y of Sharov fc Liutvi ll 19831), 
|Sh arov. Lvutvi. fc Ikonnikoval }l993~ and iBattistini et al.l 



1 19871 ). We therefore believe that our photometry of this 
cluster is reliable. 

We note that, for other clusters our p h otom etry agrees 
well with that of iReed. Harris, fc Harrij (|l992l ) and that 
B383 is brighter in all of their observations (B, V and R 
bands). We therefore believe that this discrepancy is unlikely 
to be due to an error in their photometry. This raises the 
possibility that the cluster luminosity may have genuinely 
varied between our observations. The increase in luminosity 
of B383 could have been produced by a transient in the clus- 
ter. To explain the observed variability, this transient would 
have to have bluer colours than the cluster and a brightness 
of Mv ~ -7.9. 

A potential candidate for this increase would be 
a classical nova in the c luster at the time of the 
iReed. Harris, fc Harrid l| 19921 ) observations. Novae have typ- 
ical luminosities of —6 < My < —9 and co u ld ex plain 
this blue excess in the iReed. Harris, fc Harris! (|l992l ) ob- 
servations. A classical novae of this brightness would be 
expected to have a very short outburst duration (e.g. 
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IWarnej [l995l ') and would therefore be expected to have 
faded by the time of our ob s ervati ons and even those of 
ISharov. Lvutvi. fc Ikonnikoval (jl992[ ) ~10 months later. 

Globular clusters are expected to host classical no- 
vae. There is evidence for Classical novae in the Galac- 
tic G Cs M80 l|Pogson|[l86^ IWehlau. Butterworth. fc Hogd 
Il990h and (possibly) M14 {Hogg 1964; Margon ct al. 19!^ 
Classical nov ae have also been detected in a GC in M87 
JShara et all 12004.) a nd two of M31's other GCs [Bill: 
Quimbv et aD (l2007l) : IShafter fc Quimbvl (|2007l ) and B194: 
Henze et al.l l|2009h ]. Confirmation of a classical nova in B383 
is very difficult as any remaining signatures of the event will 
be very faint. However, it offers a plausible explanation for 
such a large brightness variation. 

2.3 NIR Photometry 

2.3.1 K-band data 

To obtain K-band photometry of MSl's GCs, images across 
M31 were obtained using the Wide Field CAMera (WF- 
CAM) on the UK Infrared Telescope (UKIRT) under the 
service program USERV1652. The large field of view of the 
WFCAM makes it ideal for such a project. The coverage 
of these data is shown in figure [1] They do not currently 
cover the whole GC system, missing both the most central 
and most distant clusters. T he details of these obs ervations 
were originally presented in IPeacock" et al.l (|2009l ) but are 
summarised again below. 

The data were taken on the nights of 2005 November 30 
and 2007 August 06 with K-band seeing of 0.85-1.00" and 
0.6-0.8" respectively. To ensure the images were well sam- 
pled, each observation was taken with 2x2 microstepping to 
give an effective pixel size of 0.2". Five observations were 
taken of each field giving a total exposure time of 225s and 
a 3(T detection limit of ~19 mag. 

The images w ere reduced usin g the standard WFCAM 
pipeline (see e.g. iDve et al] |2006| ). The pipeline process- 
ing reduced and stacked the raw images and interlaced 
l|Fruchter fc Hoola 2002) the microstepped images together. 
The pipeline also applies an accurate astrometric solution 
to the images based on matching sources to the 2MASS 
catalogue. This met hod has been sh own to give positions 
accurate to 80mas (|Dve et al.l |2006| ). We determined the 
photometric zero point for each observation by calibrating 
against the 2MASS catalogue. This was done by compar- 
ing instrumental magnitudes of bright, unsaturated, stars 
in each field with the 2MASS Point Source Catalogue. This 
places the K-band photometry on the standard 2MASS 
(Vega-based) photometric system. This method has previ- 
ously been shown to give zero p oints for K-band WF CAM 
images to better than 0.02 mag (|Hodgkin et al.ll2009l) . 

2.3.2 Photometry 

Photometry was obtained for all GCs and candidates with 
WFCAM images using SExtractor. SExtractor was run in 
the same way used to obtain the optical magnitudes (de- 
scribed in section 2.2). The aperture required to determine 
the total K-band luminosity was again selected for each clus- 
ter from curves of growth. The aperture size was selected 
independently of the aperture used to determine the total 
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Figure 4. Comparison with previous K-band photometry from 
the RBC {bottom) an d profile fits to these WFCAM images from 
iPeacock et al. 1 1 I2OO9I) {top). The errors quoted are from our pho- 
tometry only, as errors are not available from the previous work. 

(/-band magnitude of the cluster. In many cases a smaller 
aperture was required to enclose all the K-band light of the 
cluster, with an average aperture size of ~4.6" selected and 
84% of the apertures 6". The use of smaller apertures for 
the K-band images is expected because of the smaller PSF 
of these images. 

As with the ugnz photometry, it was found that the 
error in the K-band photometry was often dominated by 
non statistical errors. The zero point calibration error of the 
WFCAM images is estimated to be 0.02 mag. This is larger 
than the statistical errors for most of the clusters. The K- 
band luminosity of the clusters also suffers from errors on the 
background estimation and contamination from neighbour- 
ing sources. We estimate the effect of this on the accuracy of 
our photometry using the same method used for the total g- 
band magnitude (by retaking photometry through a larger 
aperture). The median estimated error, due to contamina- 
tion, was found to be 0.025 mag. The estimated error due 
to contamination and background variation was combined 
with the statistical and calibration errors and quoted as the 
final error OK,tot in table 1. 

2.3.3 Comparison with previous photometry 

The most complete NIR data currently available for M31's 
GCs is that from the RBC. This includes K-band magni- 
tudes of 279 co nfirmed GCs i n the RBC obtained from the 
2MASS archive (|GaUeti et al. , 2004). This 2MASS photome- 
try is from either the 2MASS Point Source Catalogue or Ex- 
tended Source Catalogue and measured through apertures 
with radii of 4 and 5" respectively. The bottom panel of 
figure |4] compares our K-band photometry of all confirmed 
GCs with the K-band photometry in the RBC obtained from 
2MASS. Errors are not included for the K-band photometry 
in the RBC, so only the total errors on our photometry are 
included in these plots. 

Two clusters are found to have very different magni- 
tudes and lie off this plot. One of these (B090) is very faint 
in the previous photometry and has a very blue J-K colour. 
We therefore believe the previous photometry for this ob- 
ject is unlikely to be accurate. The other (B041) is found 
to be fainter in our photometry. The WFCAM and 2MASS 
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Figure 5. WFCAM and 2MASS images of B041. Both images 
are 45" X 45" and demonstrate the improved spatial resolution 
and signal to noise of the WFCAM images over 2MASS. 

images for this cluster are shown in figure [S] This compari- 
son demonstrates the superior depth and spatial resolution 
of the WFCAM images over 2MASS. The circles show the 
aperture sizes used for our photometry and the 2MASS pho- 
tometry. The improved spatial resolution helps to separate 
clusters from nearby sources and allows the use of signif- 
icantly smaller apertures. It is clear from this image that 
the previous photometry for these faint sources is unlikely 
to be as reliable as that presented here. From examination 
of the 2MASS catalogues we can only identify this source in 
the 'reject' catalogue. We therefore believe that the 2MASS 
photometry of this cluster is unreliable. 

The errors on the 2MASS photometry are expected to 
be larger than those obtained here since 2MASS is signif- 
icantly shallower than our data. Taking this into account, 
most of our photometry is found to be consistent, although 
there are several outliers. We believe that most of these dif- 
ferences are due to the improved spatial resolution of our 
data over 2MASS which makes it easier for us to remove con- 
tamination from nearby sources and to estimate the back- 
ground level more reliably. The poorer resolution could re- 
sult in both overestimation of the cluster magnitudes (if 
nearby sources are included in the cluster aperture) and un- 
derestimation of their magnitudes (if unresolved background 
sources result in an overestimation of the background level) . 
This highlights the importance of spatial resolution, even 
in obtaining integrated magnitudes. The brightest clusters 
are also found to be ^ 0.05 mag brighter in our photom- 
etry than in the RBC. We believe this is due to our use 
of larger apertu res for larger clusters. This was identified 
and discussed bv lGalleti et al.l l|2004l ') who attempt to apply 
aperture corrections to these clusters. However, we believe 
our use of larger apertures should be more reliable. 

An alternative method to aperture photometry is to fit 
the profile of the clusters and find their integrated mag- 
nitudes. This method removes aperture affects because it 
integrates the magnitude out to the tidal radius of the clus- 
ter. It also accounts for contamination from nearby sources 
since it assumes the cluster to have a smooth profile. This 
provides a very useful independent method of estimating the 
total magnitudes of the clusters. The results of fitting the 
profiles of clus t ers in these images are already presented in 
iPeacock et all l|2009l) . The top panel of figure |4] compares 
this integrated K-band magnitude with the aperture mag- 
nitudes found here. Again errors are not available for the 
profile fit magnitudes but are expected to be of a similar 



size to the errors obtained from aperture photometry. Three 
of these clusters are found to be brighter in our photom- 
etry. Examination of these clusters revealed that they all 
have very bright nearby neighbours. Since these cause sig- 
nificant background gradients across the cluster profiles, we 
believe that the model fits to these will be less reliable and 
the aperture photometry is probably more accurate. Some 
of the fainter clusters are also found to lie slightly outside 
2(T. As discussed later, we believe this is due to the King 
model fits being less reliable for these faint clusters. For 
these faintest clusters it is likely that aperture photometry 
gives more accurate magnitudes. 

The scatter in these comparisons highlights the diffi- 
culty in determining the NIR magnitudes of clusters pro- 
jected onto stars and surface brightness fluctuations from 
M31. We believe our approach gives the best estimate of 
their magnitudes and the most realistic errors to date. In to- 
tal we present K-band photometry for 319 and 603 sources 
classified as confirmed and candidate clusters in the RBC 
respectively. This includes the first K-band photometry for 
126 confirmed clusters and 429 candidate clusters. 



2.4 Summary of photometry 

The ugriz colours, total g and total K-band luminosity of 
M31's GCs and candidates are presented in table 1. This ta- 
ble includes the statistical errors on the ugriz colours and the 
errors on the total g and K-band luminosity (which include 
the calibration errors and estimated error due to background 
variation and contamination from nearby sources). 

It should be noted that the ugriz photometry presented 
here is on the standard SDSS (AB) photometric system, 
while the K-band photometry is on the standard 2MASS 
(Vega-based) photometric system. The magnitudes can be 
converted betw een the two systems using the following off- 
sets taken from lHewett et al] (|2006l ) : 



'iVcga = ItAB — 0.927 (5) 

gVoga = ffAB + 0.103 (6) 

rvcga = rAB -0.146 (7) 

2Voga = lAB — 0.366 (8) 

ZVcga = ZAB — 0.533 (9) 

Kab = Kvcga + 1.900 (10) 



The names of the objects in table 1 are taken from the 
RBC. The positions of the sources are taken from their loca- 
tions in our r-band images and should be accurate to better 
than 1". Some of the proposed clusters are not detected 
(or not located) in the SDSS images. The names, locations 
and classifications of these clusters (taken from the RBC 
or C09) are included in table 1. This table lists all previ- 
ously proposed clusters and candidates in the RBC. Many 
of these objects are found by this (and other) studies not 
to be genuine clusters. The classifications of these sources 
are discussed in the next section. Only those objects with 
classification flag, f=l should be considered confirmed old 
GCs. 



Table 1. Photometry and Structural parameters of M31 GCs 



GC Name^ RA^ DEC^ Classification^ Photometry Structural parameters* 
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.20 


2, 


.56 


21, 


.8 


4.11 


B028-G088 


10. 


,31871 


40, 


.98422 


1 


1 


old 


5.2 


17 


.31 


1, 


,44 


0.66 


0.38 


0.24 


0.03 


0. 


,04 





.02 





.02 


0. 


,03 


4.0 


14. 


,24 


0, 


.13 


1, 


.23 


2, 


.41 


21, 


.0 


3.68 


B029-G090 


10. 


,32433 


41 


.00640 


1 


1 


old 


6.4 


17, 


21 


1 


,98 


0.88 


0.50 


0.35 


0.04 


0. 


,06 





.02 





.02 


0. 


02 


5.8 


13. 


,24 


0, 


,06 


1 


.57 


4, 


.16 


49, 


1 


3.89 



.9 



Photometry and structural parameters of all clusters and candidates studied. Note that this includes all previous clusters and candidates - the classification flag (f) should be used to select genuine clusters 
from the table. The table is available in full as Supplementary Material in the electronic edition of the jom-nal, from the VizieR archive or from http://www.astro.soton.ac.Uk/~m.b.peacock/m31gc.html. 

^ Names taken from the RBC 

■^Position of object in SDSS r band image [J2000, degrees] 

^Classification of source (sec section 3.3 for details): 1-old globular cluster; 2-candidate cluster (21-candidate old cluster, 23-candidate young cluster); 3-young cluster; 4-galaxy; 5-HII region; 6-stellar source 

fR.BC - classification from the RBC v3.5 
fc()9 - classification from C09 
'^Aperture size used to obtain total magnitudes [arcsec] 
Error on the total magnitude, includes the statistical, calibration and systematic errors 
Structural parameters of clusters from Peacock ct al. (2009) 
'^The concentration parameter [c = log(rt/r())] 

^The half light and tidal radii [parsecs; assuming the distance of M31 to be 780kpc: McConnachie et al. (2005)] 
^The cluster core density [L^q/pc^] 



10 M. Peacock et al. 



RBC class 1 



RBC class 2 



K07 class A 



K07 class B+C 



C09 old 



1 

0.8 
0.6 
0.4 
0.2 





C09 young 



— 1 — I — r 



-1 1 1 r- — I 1 — 

+ +++++ 



. ! Am* 



— I 1 1 1 1 1 r- 



++ % 



it. 



I III ^rtd^H ji^H ii^'i^i 

_i i' '"■"TffT'llflillinihh" 



ill! "f^iiiiii^iii*^! 



■ . + t""trt - i tiil i 



11 12 13 14 15 16 17 
K [mag] 



11 12 13 14 15 15 17 
K [mag] 



11 12 13 14 15 16 17 
K [mag] 



11 12 13 14 15 16 17 

K [mag] 



11 12 13 14 15 16 17 
K [mag] 



11 12 13 14 15 16 17 

K [mag] 



Figure 6. Stellarity of objects classed as confirmed clusters (red) and candidate clusters (grey) in the RBC (their class 1 and 2 sources 
respectively) and from K07 (their class A and B /C sources respectively) . Also included are the re-classifications for many of these sources 
from C09 (their old and young clusters). The stellarity is based on SExtractor photometry of the WFCAM images. Extended sources 
have a stellarity close to and point sources close to 1. 



3 CLASSIFICATION OF SOURCES 
3.1 Stellarity 

The WFCAM images of M31 have a PSF of 0.6-0.95" corre- 
sponding to a spatial resolution of 2.1-3.6 pc at the distance 
of M31. This is a significant improvement over most of the 
images previously used to classify clusters. It allows us to 
investigate possible contamination in the previous GC cata- 
logues from single stars and previously unresolved asterisms 
of stars. Figure [6] shows the SExtractor K-band stellarity 
flag for confirmed and candidates clusters in the RBC (left) 
and, separately, the clusters and candidates from K07 (mid- 
dle) . Also included is the stellarity of old and young clusters 
from C09 (most of which are re-classifications of sources in 
the other two catalogues). 

The majority of sources can be identified as either hav- 
ing stellar profiles (with a stellarity close to 1) or extended 
profiles (with stellarity close to 0). It can be seen that some 
objects with K>15 have uncertain stellarity flags. The abil- 
ity of SExtractor to determine the stellarity of a source is 
mainly dependent on the signal to noise of the source, the 
PSF of the image and crowding around the source. From 
visual examination of the sources with uncertain stellarity 
flags, it was found that the majority of them have nearby 
sources contaminating their profiles. In general we consider 
objects with a stellarity <0.4 to be extended. However, it 
is clear that this classification is less reliable for those ob- 
jects with uncertain flags. For these objects we rely on visual 
examination of the cluster to estimate their nature (as de- 
scribed in the next section). 

It can be seen from flgure [6] that excellent agreement 
is found between our data and the classiflcations of C09 
with all sources they classify as old being extended. We also 
flnd that the majority of confirmed clusters in the RBC are 
extended. However, there are 12 RBC class 1 objects which 
have either stellar or uncertain stellarity fiags. We note that 
some of these clusters have already been reclassified by C09 
as stars. It can also be seen that many of the sources classed 
as confirmed clusters by K07 are found to be unresolved. 
We note that their work was based on images with poorer 
spatial resolution and we reclassify many of these objects as 
being stellar sources. 

Some of the young clusters from C09 are found to be 
extended and look like normal centrally concentrated GCs. 



However, it can be seen from figure [S] that many of the 
proposed young clusters have stellar, or uncertain, stellar- 
ity flags. This is likely because these young clusters can 
appear as resolved asterism s in the K-band images. This 
has pr eviously been noted bv lCohen. Matthews, fc CameronI 
(|2005l ) who used K-band images taken with adaptive optics 
to demonstrate that 4 proposed young clusters may be aster- 
isms. However, as discussed by C09, young clusters are gen- 
erally faint in K and may be dominated by only a few bright 
(resolved) supergiants making them appear as resolved as- 
terisms of stars, rather than an extended cluster. Many of 
these objects have subsequently been confirmed by HST im- 
ages to be genuine clusters. We therefore do not reclassify 
any of the proposed young clusters which appear as resolved 
stellar sources in our K-band images. 

Our data also allow us to classify many of the previ- 
ously unclassified candidate clusters. In total we classify 368 
previous candidates as likely to be stellar sources. For the 
above reasons, it is possible that we may potentially include 
some genuine young clusters in this classification. Figure [6] 
demonstrates that a large group of the proposed candidates 
are extended in our images. These objects are therefore likely 
to be either genuine clusters or background galaxies. These 
candidates represent ideal targets for foUowup spectroscopy 
in order to confirm their nature. 



3.2 Visual examination 

As discussed above, the stellarity flag for some of the fainter 
objects is relatively uncertain. Visual examination of these 
objects can help in deciding whether they are extended or 
stellar sources. Visual examination of the clusters and can- 
didates which are confirmed as extended also provides a 
method of identifying background galaxies. These were iden- 
tified as either having spiral structure, or extended elliptic- 
ity. While this method is relatively subjective, it is helpful 
in classifying an object. Inspecting the images of the objects 
also provides a useful check on our otherwise automated 
classifications. During this process, we also ensured that our 
automated photometry had selected a reasonable aperture 
size for each cluster, in order to measure its total luminosity. 
We examined the ugriz and K-band images of every 
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Figure 7. Colours of clusters and candidates from the RBC (left), K07 (middle) and C09 (right). Included are objects classified as 
confirmed old clusters (red), young clusters (blue), and candidate clusters (grey) from each catalogue. The black points show the colours 
of Milky Way GCs. The arrow represents a reddening of E(B — V) = 0.1. 



cluster and candidate studiecQ. We first examined the ob- 
jects in our sample wliicli liave recently been classified from 
tiie spectroscopic study of C09 as being background galax- 
ies. We then examined the previously classified confirmed 
clusters, followed by the proposed candidate clusters. In this 
way, we were able to reclassify some of the objects based on 
their appearance. We note that some of the newly confirmed 
galaxies look very similar to typical GCs. This highlights the 
limitations of visual examination on identifying galaxies. We 
do not reclassify any of the previously confirmed clusters as 
galaxies based on this visual examination. However, we did 
identify 3 candidate clusters with clear spiral structure and 
30 other candidates which are likely background elliptical 
galaxies (this is in addition to the candidates confirmed to 
be galaxies from the spectroscopic study of C09). During 
this visual examination it was also found that some of the 
clusters and candidates with uncertain stellarity fiags from 
our SExtractor photometry are likely to be multiple stellar 
sources, rather than extended clusters. 



^ Thumbnail images of these clusters are available at 
|http://www.astro.soton.ac. uk/~m.b. peace ck/m31gc. html I 



3.3 Colours 

Figure [7] shows the colours of objects previously classified 
as confirmed GCs (red) and candidate GCs (grey) in the 
RBC (left) and by K07 (middle). The right panel shows the 
colours of the objects which are confirmed by C09 to be 
old clusters. Shown in blue are the proposed young clusters 
from C09 and the confirmed clusters from the RBC which 
are flagged as being potentially young. For comparison, the 
black points indicate the colours of the Milky Way's GCs. 
The /?-K colour fo r the Milky Way GCs were taken from 
ICohen et al" ("2007^ and optical colours from the Harris cata- 
logue (Harris 1996). The colours of the Milky Way GCs were 
tr ansformed in to the ugriz filters using the transformations 
of I Jester et al.l (2005 ) and dereddened using the values for 
E{B — V) quoted in the Harris catalogue. Only the Milky 
Way clusters with E(B — V)< 0.4 are included. This limits 
the Milky Way sample to mainly low metallicity clusters. 
It can be seen that the Milky Way's GCs define a tight re- 
gion in the colour-colour plots. For this reason, the colours 
of the proposed GCs and candidates in M31 are very useful 
in classifying the objects. 

It should be noted that the colours of M31's GCs are 
reddened due to both Galactic extinction and extinction in- 
trinsic to M31. The Galactic reddening in the direction of 
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Figure 8. Top: Colours of previously confirmed clusters which are confirmed here to be extended (red). Bottom: Colours of proposed 
candidate clusters which are confirmed to be extended (grey). The lines indicate linear fits to the colours of all confirmed M31 GCs. The 
black points indicate the colours of the Milky Way GCs and the arrow represents an extinction of E{B — V) = 0.1. 



M31 is relatively uncertain, but it is estimated for the re- 
Rion around the disk of M31 to he E(B - V) ~0.062 mag 
ijSchlegel. Finkbeiner. fc Davij lTogsl). However, the extinc- 
tion due to M31 itself can be much larger and varies sig- 
nificantly between GCs due to their locations in (and line 
of sight depths through) the galax y. Previous work (e.g. 
iBarmbv etallbOOOl : iFan et al .1120081 ) has demonstrated that 
the reddening for some of these clusters can be substan- 
tial. For example, the very red cluster in figure [7] with g- 
K=6.95 is B037 which is known to b e heavily reddened 
[E{B = 1.38: iBarmbv et all (|2000l )]. 

Figure [7] shows that our colours are in good agreement 
with the classifications of C09. It can be seen that most of 
the objects classified by C09 as old clusters define a tight re- 
gion which is consistent with the (reddened) colours of the 
Milky Way's GCs. In most cases the confirmed clusters in 
the RBC also have colours consistent with the Milky Way's 
GCs. Many of the confirmed clusters from K07, and a few 
of the confirmed clusters from the RBC, have colours which 
are not consistent with the Milky Way's GCs or the ma- 
jority of the confirmed GCs in M31. This is in agreement 
with our conclusions from the previous section that some of 
the previously confirmed clusters may be misclassified stars. 



The colours also suggest that many of the unclassified candi- 
date clusters may be stars, asterisms of stars or background 
galaxies. 

3.3.1 Young clusters 

These colours clearly identify the population of very blue 
clusters that have been noted by previous studies. It can 
be seen that our colours are in excellent agreement with 
the spectroscopic classifications of C09. We also find good 
agreement with the confirmed clusters in the RBC which 
are fiagged as potential young clusters (fla g yy=l,2 or 3 in 
the RBC). This flag is based on the work of lFusi Pecci et ahl 
(2005). Most of the previously identified young clusters are 
much bluer in g-r t han any GC in t h e Mil ky Way. Using 
a similar method to iFusi Pecci et al.l (|2005l ). we define all 
objects with i;-r<0.3 to be young clusters. 

Some of the proposed young clusters have colours which 
are consistent with being old clusters. However, these objects 
are also consistent with being young clusters with reddened 
colours. These clusters are also found in high density regions 
of M31 and look similar to the other young clusters we have 
observed. We therefore choose to keep the previous (spectro- 
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Table 3. Classifications of sources 



Classification 


Number in 
this study 


RBC 1 


Previous classification of these objects 
RBC 2 K07 A K07 B/C C09 old 


C09 young 


1: old globular cluster 


416 


342 


41 


27 





336 





2: candidate cluster 


373 


6 


101 


9 


256 


3 





3: young cluster 


156 


46* 


78 


2 





2 


151 


4: background galaxy 


189 


5 


170 


4 


10 








5: HII region 


17 





14 


3 











6: stellar source 


444 


10 


153 


66 


215 


1 





Total (previous 


catalogues): 


409 


557 


111 


481 


342 


151 



* Many of these RBC class 1 clusters are flagged separately in the RBC as potentially young clusters. 



scopic) classification for these clusters and suspect that their 
colours may be reddened. It can also be seen from figure [7] 
that two clusters classified as old by C09 (B386 and PHF7- 
1) have very blue colours. We reclassify these two objects as 
young clusters. 

3.3.2 Old globular clusters 

Figure [8] shows the colours of all confirmed and candidate 
clusters following the removal of all stellar objects based 
on their stellarity flag or visual examination of the cluster 
images. We have also removed those objects identified in 
the previous section as being as young clusters. It can be 
seen that, having removed these objects, the colours of the 
confirmed clusters are now consistent with the colours of 
the old GC system of the Milky Way. The clusters extend to 
much redder colours, but this is consistent with the expected 
reddening due to extinction from M31. 

The grey points in the bottom panels of figure |8] show 
the colours of the remaining candidate clusters after the re- 
moval of non-extended objects. Comparison with the con- 
firmed old clusters shows that many of the candidates have 
colours consistent with being old clusters. These clusters are 
fiagged as old candidates in table 1 and should be considered 
the strongest candidate clusters. We also identify candidates 
with very blue colours, consistent with the other young clus- 
ters identified. These are flagged as young candidate clusters 
in table 1. It can be seen that, despite removing objects iden- 
tified as stars, the colours of many of the candidate clusters 
are inconsistent with being either old or young clusters. As 
we are uncertain of the classification of these objects, we 
retain their classification as candidates. However, it is likely 
that many of these candidates are either background galax- 
ies or unresolved asterisms. 



3.3.3 Extended clusters in the halo of M31 

Recent studies of the halo of M31 have identified a 
population of very extended cluster s ijHuxor et al ] I2OO5I : 
iMackev et alll200g : llfuxor et al.ll2008h . These clusters have 
half light radii much greater than the majority of clusters in 
M 31. For a de scription of these clusters we refer the reader 
to iHuxor et a l. (2008). Seven of these clusters are located 
in our SDSS images and can be identified in table 1 from 
their names which have the prefix HEC ('Halo Extended 
Cluster'). Our colours of these clusters were found to be less 
reliable than the other clusters studied. This is because they 



are resolved, due to the diffuse nature of the clusters, into 
multiple sources. 

The colours for these clusters were therefore re- 
measured through 12" apertures using the IRAF:APPHOT 
task PHOT. A smaller aperture of 8" was used for HECll 
due to a bright neighbouring star. This method gives reli- 
able results for clusters in the halo of M31 where there is lit- 
tle contamination from neighbouring sources and the back- 
ground is relatively smooth. None of these extended clusters 
are identified in the inner regions of M31, although detecting 
such extended and faint objects in front of the M31 would 
be very difficult. 

These clusters are identified in figure |8] as open green 
points. It can be seen that the colours of these clusters are 
now consistent with the other old GCs in M31. The errors on 
the colours of the HEC clusters are larger than those of the 
other GCs. This is due to their diffuse nature and the use of 
large apertures, which increases the total sky background. 

3.4 Final classification 

Our final classification is based on: the stellarity of the 
object; its colours; visual examination of the object in our 6 
bands; velocity information and classifications from previ- 
ous studies. Table 1 lists these classifications for all GCs and 
candidates. For comparison we also include the previous 
classifications from the RBC and C09. For consistency we 
have tried to keep our classifications similar to those used 
in the RBC. If we have no reason to reclassify the sources, 
we keep the original classifications (where available from 
C09, which were found to agree best with our classifica- 
tions, otherwise from the RBC). The classifications used are: 

1: old globular cluster: extended and has colours con- 
sistent with the Milky Way's GCs. Its velocity is confirmed 
from previous work (K07, C09 or the RBC) to be consistent 
with being in the M31 GC system, or the object is confirmed 
from high resolution HST images. 

2: candidate cluster: not confirmed, but previously proposed 
as being a cluster or candidate and is found here to be 
extended (or have uncertain stellarity). Candidate is sub- 
divided, depending on whether its colours are consistent 
with being an old cluster (21), consistent with being a 
young cluster (23) or inconsistent with being a cluster (2). 
3: young cluster: has colours consistent with being young. 
If previously classified, may appear as a resolved asterism 
in K, but looks like a cluster in the SDSS images. 
4: background galaxy: previously classified from spec- 
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Figure 9. K-band GCLF for all sources classed as confirmed GCs 
(solid, red) and young clusters (open, blue). 



troscopy by C09 or identified from our visual examination. 
5: HII region: from previous classification of C09. 
6: stellar source: object appears to be a single stellar source 
or a previously unresolved asterism of stellar sources. 

The total number of sources of each class are shown in 
table |3l For reference we include whether these objects were 
previously classified as: clusters or candidates in the RBC 
(RBC 1 and RBC 2 respectively); clusters or candidates by 
K07 (K07 A and K07 B/C respectively); old or young clus- 
ters by C09 (C09 old and C09 young respectively). It can be 
seen that we have reclassified 10 previously confirmed clus- 
ters in the RBC as likely stellar sources. We also reclassify 
6 of these objects as candidate clusters, as we are uncer- 
tain of they nature, or they lack spectroscopic confirmation. 
Some of the candidate clusters in the RBC are confirmed 
to be old or young clusters. This is based on the new spec- 
troscopic confirmations by C09. We are also able to classify 
many of the candidate clusters in the RBC as stars. In most 
cases we find good agreement with the new classifications 
of C09. Their catalogue includes fewer objects because they 
do not provide classifications for the whole GC system. All 
objects classed as young clusters by C09 are retained in our 
classification. 

We have reclassified many of the confirmed clusters 
from K07 as likely stellar sources. We have also been able 
to classify nearly half of the cluster candidates from this 
catalogue as being stellar. We believe this is due to our im- 
proved spatial resolution compared with the images used for 
this previous catalogue. We identify and remove 8 objects 
from the catalogue of K07 which are within 2" of another 
previously identified object in the RBC and we believe are 
now duplicated in the RBC. A further 5 objects from the 
catalogue of K07 appear to be associated with objects in 
the catalogue of C09. The names for these objects in table 
1 are the combination of their identifications in each cata- 
logue. 



4 PROPERTIES OF CONFIRMED GCS 

Figure [9] shows the GC Luminosity Function (GCLF) for 
all confirmed GCs (solid bars) and young clusters (open 



bars) with K-band photometry. These clusters are not cor- 
rected for extinction. However, extinction is not very sig- 
nificant in the K-band where the maximum correction for 
the most extreme case of B037 is only 0.5 mag (the width 
of the bins used). The peak of the GCLF is found to be at 
K~14.2 mag. The K-band luminosity of a cluster is a useful 
estimate of its mass. This is because, in addition to being 
less effected by extinction, the K-band mass to light ratio 
{M/L) is less effected by metallicity than optical bands. The 
mass to light ratio of a 12 Gyr cluster in the K-band has 
previously been estimated to be 0.9<M/Z/<1.3 for metal- 
licities in the rang e 0>[Fe/H]>-2 (|Bruzual fc Charlot|[200l : 
iForbes et al. I l2008i ). To estimate the peak mass of the old 
GCs in M31, we assume a K-band M/L ratio of 1.1 for all 
clusters (as the metallicities are not kno wn for all of the clus- 
ters). At the distance of M31 (780 kpc: iMcConnachie et al.l 
2oog) and assuming the K-band magnitude of the sun to 
be MksP)=3.29 rnag [Mkr=3.33: |Coj (l2000|); Mk.0=Mk0- 
0.04: ICarpenteil (|2003)], this implies a peak mass of Mpeak ^ 
3 X IO^Mq. Th is is sUghtly higher than that found for Milky 
Way GCs (e.g. ICohen et al.|[2007l ). However, this difference 
is relatively small compared with the expected uncertainty 
in the peak mass. This is due to errors in accurately esti- 
mating the peak in the GCLF combined with errors on the 
distance to M31 and the value used for the mass to light 
ratio. 

For the fainter GCs, it is likely that masses esti- 
mated from their integrated K-band luminosities are less 
accurate due to stochastic effects. Stars at the tip of 
the red giant branch at the distanc e of M31 are ex- 
pcctcd to reach mag nitudes of K=17.5 (|Ferraro et al.|[200Gl : 
pTabur, Kiss, fc Bedding 2009). While stars this bright are 
relatively rare, the integrated light of some fraction of these 
faint clusters can therefore be dominated by a relatively low 
number of these stars. It is also likely that some of the 
faintest clusters in M31 are missing from our catalogue. 
These clusters should be detected in our data. However, 
identifying these faint clusters in front of M31 would be very 
difficult. 

As expected the proposed young clusters peak at fainter 
magnitudes than the old GCs. Some of these clusters are 
found to be relatively bright, reaching luminosities similar 
to the peak of the GCLF. This suggests they are more mas- 
sive than typical young open clusters in the Milky Way. 
Whi le this is in agreement w ith the conclusions of other work 
(e.g. lFusi Feed et ahllioosl ). it should be noted that our con- 
clusions based on this K-band luminosity are limited. The 
M/L ratio of these clusters is likely to be significantly lower 
than the M/L ratio for the older clusters in the galaxy. Also 
stochastic effects in these young and faint clusters are likely 
to be significant in the K-band. 



4.1 The structure of M31 GCs 

The excellent spatial resolution of these WFCAM images 
allows us to investigate the structure of M3rs old GCs. This 
can be done by fitting PSF convolved King models to their 
profiles. The results of this fitting to 239 class 1 c l usters 
from the RBC are already presented in lPeacock et al.l (|2009l ) 
(hereafter P09). For convenience, we have included these 
structural parameters, where available, in table 1. Here we 
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Figure 10. Structural parameters of M31 GCs {top) and Milky Way GCs {bottom). The dashed line indicates the GCs in the Milky 
Way over a similar galactocentric radius. 



review these data and use it to consider the structure of 
M31's GC system. 

From our new classifications it is found that 186 of the 
239 clusters studied by P09 are confirmed as old GCs. We 
have also determined the structural parameters for an ad- 
ditional 27, newly confirmed, GCs with WFCAM images. 
These parameters were determined using exactly the same 
method as P09. This results in a sample of the 213 old GCs 
with structural parameters (representing over half the con- 
firmed GCs identified in M31). The reliability of these struc- 
tural parameters, was investigated by P09 and they were 
found to be consistent with those obtained for 33 clusters 
whose parameters had already been measured using higher 
spatial resolution HST images (iBarmbv et al.1 2007). They 
also found that the structural parameters obtained from fit- 
ting different WFCAM images of the same GC were consis- 
tent. However, the reliability of the parameters was found 
to decrease significantly for clusters with K>15 mag. We 
note that this decrease in reliability for the fainter clusters 
is in agreement with the si gnal to noise limit proposed by 
ICarlson fc HoltzmanI l|200ll '). The half light radii of the clus- 
ters was found to be the most reliable paramet er (as found 
by other studies e.g. iKundu fc Whitmorelll998l ). Errors are 
not available for these parameters, but we believe them to 
be relatively reliable for most clusters with K<15 mag. 

Figure [To] shows the concentration [c =log(rt/ro)], core 
radii (ro), half hght radii (r^), tidal radii (rt) and core den- 
sity (po) of M31's old GCs. For comparison the same pa- 
ram eters for the Milky Way's GCs are shown on the bottom 
row (|Harrislll996l ). Before comparing these populations there 
are two important differences which need to be considered. 
Firstly, the M31 parameters are based on the K, rather than 
the V-band luminosity of the clusters. As a result we expect 
offsets in parameters such as the core luminosity density 
due to the different mass to hght ratios. However, this is 
unli kely to have a significant effect on the size of the clus- 
ters l|Cohen et a l.l l2007i ') . Secondly, due to a lack of WFCAM 
data, the sample of M31 GCs does not include the most cen- 
tral or distant GCs in the galaxy. It can be seen from figure 
[To] that there are a lack of core collapsed GCs in our sample 
compared with the Milky Way. This can be partially ex- 



plained by the exclusion of the innermost GCs (where most 
core collapsed GCs are located in the Milky Way). Taking 
this into account, we still find fewer of these clusters than 
expected. Potentially some faint core collapsed clusters may 
be missed by GC surveys as they would be the most difficult 
clusters to resolve. However, it is likely that these clusters 
are present, but have their concentrations underestimated. 
This is because their core radii will be much smaller than 
the PSF of our images, making it very hard to deconvolve 
and measure them. Comparison of the core radii of the clus- 
ters does show that we are missing, or overestimating, the 
core radii of some of the very smallest core radii clusters. 
A lack of core collapsed clusters can also be s een in similar 
profi le fits to HST images of Cen A clusters (|Jordan et al.l 
|2003). Allowing for these effects, we see no strong evidence 
for differences between the structure of the old GCs in M31 
and Milky Way. 

In the Milky Way, it is known that some GC properties 
are related to position in the Galaxy. In figure [TT] we plot 
the structural parameters of M31's GCs as a function of 
their projected galactocentric radius {Rgcp '■ taken from the 
RBC). To help identify potential trends in these data, we 
have binned the clusters into groups of 25 and determined 
the median value of their parameters for each group (bold 
points). It can be seen that the luminosity of the GCs ap- 
pears to decrease slightly with Rgcp- This is mainly driven 
by a deficit of faint clusters in this region combined with the 
most massive clusters being centrally located. The relatively 
low number of faint clusters in the central regions is likely 
due to selection effects, as it is very difficult to identify these 
clusters projected against the dense central regions of M31. 
It is therefore likely that some of these clusters are miss- 
ing from our catalogue. Selection effects can not explain the 
lack of very massive clusters in the outer regions. However, 
there are very few of these very bright clusters. The distri- 
bution also suggests that more central clusters have smaller, 
more concentrated cores. This relationship may be expected 
from the evolution of the GC system as central GCs are 
expected to evolve more quickly due to greater interactions 
with their h ost galaxy. This effect is ob served in the Milky 
Way's GCs l|Diorgovski fc Mevlanl[l994l ) . However, we again 
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Figure 11. Properties of M31's GCs (crosses) as a function of 
projected distance to the center of M31 (Rgcp)- Bold points show 
the median values for clusters bin ned on Rgcv The lin e demon- 
strates the relationship found by iBarmbv et al. I ll2007l) between 
Tfi and -Rgcp. 



colours are consistent with their spectroscopic classification 
by C09 as young clusters. We note that many of these clus- 
ters look like resolved asterisms in our K-band images. How- 
ever, some of these are confirmed by HST images to be gen- 
uine clusters. Higher spatial resolution optical images than 
available here are required in order to confirm their nature 
as genuine young clusters. 

We have identified that many of the confirmed clusters 
from iKim et al] l|2007h are likely stellar sources (we retain 
only 27 of their 111 confirmed clusters as old clusters). We 
also identify 10 confirmed clusters in the RBC as likely stel- 
lar sources. While we have considered the classifications from 
K07 and the RBC separately in this paper, we caution that 
all of the objects confirmed by K07 to be clusters are in- 
cluded as confirmed clusters in the current version of the 
RBC (v3.5). We also provide new classifications for many of 
the cluster candidates proposed by this previous work. We 
identify many of these candidates to be stars and reduce the 
number of unclassified candidate clusters to 357. 

Taking extinction and selection effects into account, we 
find both the colours and structure of the old M31 GC sys- 
tem are consistent with the Milky Way's. We note a po- 
tential lack of both core collapsed and very extended GCs 
in our M31 sample. We caution that some (or all) of this 
effect may be due to selection effects in identifying these 
clusters, or difficulties in accurately measuring their param- 
eters, rather than an intrinsic difference in the populations. 



caution that this is a weak trend, and that the same selec- 
tion effects could potentially prejudice us against extended, 
low density clusters in the inner regions. 

The half hght radius of the clusters can be seen to in- 
crease with Rgcp- This has previously been observed for a 
smaller number of GC s in M 31 but over a greater range of 
Rgcp bv lBarmbv et al] l|2007l). The line incl u ded in this plot 
is the relationship found bv lBarmbv et all |2003) not a fit 
to our data. This demonstrates the excellent agreement be- 
tween the trend they identify and that found here for a larger 
number of clusters . A similar trend is also found for GCs 
in the Milky Way llva,n den Bergh. Morbev. fc Pazdej|l99ll : 



iDiorgovski fc MevlanI 1994 ) and in Virgo cluster galaxies 
I Jordan et al.ll2005l ). Unlike other cluster sizes, the half hght 
radius of a cluster is thought to be largely unaffected by 
evolution. Therefore this relationship may be related to the 
properties of the globular cluster system at the time of for- 
mation. 



5 CONCLUSIONS 

Our final catalogue includes 416 old GCs. Where detected, 
we provide self consistent ugriz and K-band photometry for 
the proposed clusters and candidate clusters. We note the 
difficulty in providing accurate photometry for some of these 
clusters due to the complex background of M31. We high- 
light the need for good spatial resolution in order to remove 
contamination from non cluster light when obtaining inte- 
grated magnitudes. Where available, we find our photome- 
try to be consistent with that previously published. From 
our multicolour photometry, we confirm the population of 
very blue clusters identified previously. We show that these 
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